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ABSTRACT 



We report on the first Suzaku observation of IGR J 1 63 1 8-4848, the most extreme example of a new group of highly absorbed X-ray 
binaries that have recently been discovered by the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The Suzaku 
observation was carried out between 2006 August 14 and 17, with a net exposure time of 97 ks. 

The average X-ray spectrum of the source can be well described (^JL = 0.99) with a continuum model typical for neutron stars i.e., 
a strongly absorbed power law continuum with a photon index of 0.676(42) and an exponential cutoff at 20.5(6) keV. The absorbing 
column is Nu = 1.95(3) x 10 24 cirr 2 . Consistent with earlier work, strong fluorescent emission lines of Fe Ka, Fe K/?, and Ni Ko- are 
observed. Despite the large Nu, no Compton shoulder is seen in the lines, arguing for a non-spherical and inhomogeneous absorber. 



Seen at an average 5-60 keV absorbed flux of 3.4 x 
hours. 



10 erg cm s , the source exhibits significant variability on timescales of 
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1. Introduction 

IGR J16318-4848 was detected on 2003 Jan. 29 during a 
scan of the Galactic plane by the IBIS/ISGRI soft gamma- 
ray detector o nboard the International Gamma Ray Laboratory 
(INTEGRAL : ICourvoisier et alll2003llWalter et alll2003l) . The 
source was the first and most extreme example of a number 
of highly absorbed Galactic X-ray binaries discovered with 
INTEGRAL. Due to the strong absorption, which can exceed 
an equivalent hydrogen column of 10 24 cm~ 2 , these sources are 
extremely faint in t he soft X-rays and had not been d e tected 
by earlier miss ions dRodriguez et all l2003t iPatel et all 12004 
lKuulkersLl2005h . 

Rig ht after its discovery, a re-analysis of archival ASCA 
data by Murakami et al. (2003) revealed a highly photoabsorbed 
source (Nu = 4 x 10 23 crrT 2 ) coincident with the position 
given by INTEGRAL. The data also suggested an iron emission 
line at 6.4 keV. Thes e results were confirmed by various sub- 
sequent studies (e.g. ISchartel et all 120031: |d e Plaa et al., 2003; 
Revn ivtsev etall 120031: IWalter et all \20Q3 . iMatt & Guainazzil 
d2003l) detected intense Fe Ko-, Fe KJ3, and Ni Ko- emission 
lines in the spectrum. Based on the interstellar absorption to- 
ward the system, which is two order s of magnitude lower tha n 
the measured N n , IRevnivtsevI d2003l) . IFilliatre & Chatvl (120041) . 
and iLutovinov et all d2005l) also suggested that much of the X- 
ray absorption is intrinsic to the comp act object. 

In an optical study of the system, IFilliatre & Chatvl d2004l) 
proposed that IGR J 163 18 -4848 is a High Mass X-ray Binary 
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(HMXB) with an sgB[e] star as the mass donor surro unded by a 
dense and absorbing circ umstellar material (see also lRevnivtsevi 
120031: iMoon et alll2007l) . This dense stellar wind results in sig- 
nificant photoa bsorption within the bin ary system. Based on 
the optical data, Filliatre lfe Chatvl (120041) suggest a distance be- 
tween 0.9 and 6.2 kpc for the system. A likely locat i on for 
the source is in t he Norma-Cygnus arm (IRevnivtsevI [2003; 
IWalter et a l.. 2004), which would place it at a distance of 4.8 kpc 
(IFilliatre & Chatvl 120041) . 

In this Paper, we describe the results of follow-up observa- 
tions of IGR J 163 18 -4848 obtained with the Suzaku satellite, 
the instruments on which are uniquely suited to study Compton- 
thick absorption. In Sect. [2] we describe the data reduction. 
Sectionf3]is devoted to a presentation of the results of the spectral 
and temporal analysis. We discuss our results in Sect. [4] 



2. Data analysis 

We observed IGR J16318-4848 with Suzaku from 2006 
August 14 until 2006 August 17 for a total net exposure of 
97 ks (Suzaku sequence number 401094010). We used the stan- 
dard procedures to r educe the data from the X-Ray Imaging 
Spectrometer (XIS, Koyamaetal., 2007) and the Hard X- 
Ray Detector (HXD, Tak ahashi et all 120071) . For the XIS in 
particular we barycentered the data with aebarycen (version 
2008-03-03) and then extracted source events, images, spec- 
tra, and lightcurves with XSELECT v2.4. A circular source 
extraction region of 3'.23 radius was applied. The background 
spectrum was extracted from a circular region having the same 
area as the source extraction region. This process was done 
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Fig. 1. Spectrum of IGR J 163 18 -4848 in the range 0.3-60 keV. 



for every XIS. Response matrices and ancillary response files 
were generated using XISRMFGEN (version 2009-02-28) and 
XISSIMARFGEN (version 2009-02-28), taking into account 
the hydrocarbon co ntamination on the optical blocking filter 
dlshisaki et all 120071) . As recommended by the Suzaku team, 
the spectra of the three front illuminated CCDs (XIS0, XIS2, 
and XIS3) were then combined with addascaspec (version 1.30). 
Although the XIS1 was operational when the observation was 
made, it is not used in the present study due to cross calibration 
issues. 

To extract the HXD PIN spectrum, we again followed the 
standard procedure of barycentric correction, gti-filtered spec- 
trum extraction with XSELECT and dead-time correction with 
HXDDTCOR (version 1.50). The cosmic background was cre- 
ated with a model provided by the Suzaku team using a flat 
response (ae_hxd_pinflate2_20080129.rsp) and then combined 
with the internal background model provided by the Suzaku 
team (ae401094010_hxd_pinbgd.evt). The resulting combination 
is used for the background subtraction. The response matrix used 
for the analysis is the one proposed by the Suzaku team for the 
time of our observation, ae_hxd_pinxinome2_20080129.rsp. The 
count rates of IGR J16318-4848 are 0.1437 + 0.001 cts s" 1 for 
the combined XISs and 0.6108 + 0.004 cts s 1 for the HXD PIN 
diodes. 

For the analysis with XSPEC (v. 11.3.2ag: lArnaudLll996l) we 
rebinned the spectrum to a minimum of 250 and 200 counts per 
bin for the XIS and the PIN, respectively. The uncertainties for 
all fits are quoted at the 90% level for a single parameter of in- 
terest. In order to account for flux cross calibration issues among 
the instruments, in all spectral fits a multiplicative constant was 
introduced. 



3. Suzaku observation of IGR J1 631 8-4848 

3.1. Spectral analysis 

Although we detected a soft excess in the spectrum below 5 
keV (Fig.Q]), we did not include it in the modeling because it is 
most probably due to a se rendipitous s ource at a distance 30" 
from IGR J 163 18 -4848 dlbarra et all 120071: iMatt & Guainazzil 
120031) . The presence of this source could not be confirmed here 
because of the lower angular resolution of the XISs compared to 




Fig. 2. Confidence contours (68, 90, and 99 per cent) of the col- 
umn density and the folding energy as a function of the photon 
index. The cross mark indicates the best fit value. 




Energy [keV] 

Fig.3. Broadband spectrum of IGR J16318-4848 together with 
the best fit model and its residuals. 



XMM-Newton, even when using an optimal attitude solution for 
Suzaku by measuring the attitude directly through following the 
location of IGR J 1 63 1 8-4848 on the XIS chips. 

In order to describe the 5-60 keV broad-band spectrum of 
the source we fit the spectral continuum with an absorbed cut- 
off powerlaw, taking also into account non-relativistic Compton 
scattering. Photoabs orption was modeled with a revised version 
of the TBabs model dWilms et all l2000l 120061). using the inter - 
stellar medium abundances summarized by Wil ms et ail (|2000). 
This model describes the continuum extremely well (Fig. [3). 
In addition to the continuum, strong fluorescent emission lines 
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Table 1. Best fit parameters obtained from modeling the joint XIS and HXD data in the 5-60 keV band. 1 
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Fig. 4. Close-up of the Fe Ka band. 
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Fig. 5. Top: Lightcurve for the XIS (5-12 keV, squares) and the 
HXD PIN (12-60 keV, triangles). Bottom: Hardness ratio as a 
function of time. 



from iron (Fe Ka and K/3) and nickel (Ni Ka) are introduced in 
the model (within the absorber) to obtain a satisfactory descrip- 
tion of the data (Fig. 0). We model these lines with Gaussians 
fixed to a width of cr = 0.1 eV (i.e., we use lines narrow com- 
pared to the resolution of the XIS). The Fe Ka line is modeled 
as the superposition of the Fe Kai and Fe Ka^ lines, with the 
relative line normalizations held at the 2:l-ratio of the floures- 
cence yields of these lines and the Fe Ka2 line constrained to 
be 13.2 eV below the Fe Kai line. We also modeled the Fe K/3 
line as the combination of the Fe K/3\ and Fe K/3?, lines (the Fe 
K/3? energy being fixed to 16 eV below Fe K/3\, and its inten- 
sity to half the one of Fe K/3]). This physically correct approach 
is to be preferred to modeling the Fe Ka and Fe K/3 lines with 
a single Gaussian. We introduced a multiplicative constant c to 
normalize the HXD flux with respect to the XIS one. 

The resulting model (Table[T|) provides a good description of 
the data Cr 2 /dof = 242.6/245). With N H = 1.95+Jj ^ x 10 24 cm -2 
the column density is very high, as is to be expected for this kind 
of so urce, and is in agreement with the previous observations 
(e.g iLutovinov et all 120051: IWalter et all l2006t llbarra et all 
2007). In contrast, the photon index, T = 0.676 + ° °°^, is consid- 
erably harder th an in several earlier analyses (e.g.. IWalter et alj 
l2004t F = 2.6 or llbarra et all d2007l) : T = 1.35-1.46). As shown 
by the contour plots in Fig. [2] our broad-band data allow us to 
determine T to a high precision. The photon index is not corre- 
lated with /Vh, and there is only a slight dependency between F 
and Ef Q \d, which is much smaller than the difference between the 
photon index found here and that found in earlier observations. 



Despite the large /Vh, which corresponds to a moderately 
high Thomson optical depth of r es = 1.3, no Compton shoul- 
der is apparent in the spectrum and all lines are well modeled 
with narrow Gaussians (Fig. HJ. In order to determine an up- 
per limit_^o£_fhe_^iK inj_putative Compton shoulder, follow- 
ing iMatt & Guainazzil (|2003) we model this feature by adding 
a moderately broad (cr = 50 eV) Gaussian at 6.3 keV to the 
model. The 90% upper limit for the flux in the Compton shoul- 
der is 1.8xl0~ 5 phcirT 2 s _1 , corresponding to a 90% upper limit 
of 34.6 eV for the equivalent width. 

Data from the three XIS and the HXD-PIN were used to ob- 
tain lightcurves in the 5-12 keV and in the 12-60 keV band. To 
study the evolution of the spectral hardness of the source, count 
rates were determined at the resolution of the good time intervals 
of the XIS0 detector, which cover approximately one Suzaku- 
orbit each (~90 minutes). Figure [5] shows the significant vari- 
ability of IGR J 163 18 -4848 on this resolution. Throughout the 
observation, for XIS count rates above 0.1 counts s _1 the source 
shows no clear dependence of the hardness ratio from the source 
count rate, indicating that only slight changes in the spectral 
shape occur. At even lower count rates, the X-ray spectrum soft- 



1 We list the photon index (O, folding energy (E^m), hydrogen 
equivalent column (N H ), Fe abundance (A Fc ), the total absorbed and un- 
absorbed fluxes, and the energy (E) and equivalent width (EW) of the 
fluorescence lines. The norm of the absorbed cutoff powerlaw (A culo (f p i) 
is defined as the photon flux at 1 keV; for the absorbed Gaussian lines 
the norm (F) equals the total line flux. 
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ens, but the signal to noise in the X-ray spectrum is too low to 
allow us to quantify these changes further. 

4. Summary and Conclusions 

We have presented first results from the analysis of a long 
Suzaku observation of IGR J16318— 4848, the most extreme of 
the strongly absorbed "INTEGRAL -sources". As found in pre- 
vious studies, the average spectrum of the source is consistent 
with a strongly absorbed exponentially cutoff power-law and 
strong flourescent line emission. In contrast to earlier studies, 
the power-law photon index was found to be considerably harder 
than before (AF from 0.67 up to 1.93). This result can be due to 
the significantly better signal to noise ratio in the energy band 
above lOkeV compared to the earlier studies, which allows for 
a better determination of the high energy cutoff, the continuum 
parameters, and Ah than the earlier soft X-ray measurements, 
although an instrinsic change in the source is not ruled out. 

The soft excess below 2keV is pro bably due to a sere ndip- 
itous source near IGR J16318-4848 ( iFbarra et all 120071) . The 
considerable variability of the source can be explained as being 
due to variations in A^h. 

As pointed out by Walt er et"ai1 (|2004), the general spectral 
characteristics derived from the fit are typical for acc reting neu- 
tron stars (e.g.. iNaik & Paull 120041: Iffill et all 120081) . Note that 
this result does not mean that the neutron star nature of the com- 
pact object in IGR J 163 18 -4848 is confirmed, which would re- 
quire e.g. the detection of pulsations. A search for pulsations in 
the range between 1 s and 1 ksec was negative, while shorter pe- 
riod pulsations are probably not dete ctable due to the sm earing 
of pulsations by Compton scattering (Kuste r~et alll2005h . 

Turning to the emission lines, we note that our fit requires 
a slight overabundanc e of iron with respect to the ISM values 
of IWilms et all d2000l) . as one would expect for an evolved star. 
Furthermore, the flux ratio of Fe and Ni also points towards a Ni 
overabundance by a factor of ~2.5 with respect to Fe. 

The ratio of the Fe Ka and Fe K/3 line fluxes is given by 77 = 
(F(Fe K/3,) + f (Fe K/3 3 ))/(F (Fe K«i) + F(Fe Kor 2 )) = 0.086 ± 
0.008. This flux ratio is formally slightly smaller than that found 
in theoretical calculations for neutral gas phase Fe atoms of 
Uacobs & Rozsnvail (1 19861 n = 0.12 1). [Kaastra & Mewd d!993l 



77 0.125), or]jankowski & Polasik ( 1989, 77 0.132(2)), and it 



is also smaller than the value of 77 found in experimental mea- 
sur ements perform ed i n solid Fe (e.g., 77 - 0.1307(7) found 
bv iRai etal.lll998l and iPawlowski et all l2002h . The difference 
between the different theoretical calculations is due to certain 
approximations made in solving the structure of the excited Fe 
ion after the K-shell photoabsorption, while for the latter mea- 
surements 77 is affected by internal absorption in the Fe crystal 
used to make the measurements as well as by the dependence 
of the emission probability of the photoelectron on orientation. 
The systematic uncertainty of 77 in theory and measurements is 
therefore probably as large as 0.02, which would make our mea- 
surement consistent with neutral Fe. We note that our value for 
77 is significantly smaller than the 77 = 0.20 + °^ found in the 
XMM- Newton EPIC-pn a nalysis of iMatt & Guainazzil (I2003L 
but see I Walter et al.ll2003l) . These authors speculated that this 
higher 77 could be due to the absorbing wind being moderately 
ionized. Given that the line ratio (and also the line energy) found 
in the higher resolution Suzaku data are consistent with neutral 
Fe, we might be seeing a change in the ionization structure of 
the wind between the XMM-Newton and the Suzaku observa- 
tions. Alternatively, the larger value for 77 may be due to sys- 
tematic effects in the XMM-Newton analysis: With Suzaku, the 



Fe K/3 line and the Fe K edge are easier to separate and the 
spectral continuum is better constrained in the present analysis 
than with XMM-Newton, since spectral information is available 
above 9 keV. 

Finally, despite the large column of the source, no signifi- 
cant evidence for the presence of a Compton shoulder is found 
in the Suzaku spectrum, which is consistent with previous re- 
sults. This result is in contrast to the expectatio n for absorp- 
tion in an homogeneous medium: As shown by iMattl (|2002), 
with this assumption the equivalent width of the Fe Ka line at 
the A H of IGR J16318-484 should be much less than that ob- 
served here, and a strong Compton shoulder shou ld be present, in 
line e. g. with the Compton shoulder observed by Watanabe et all 
(20031) in GX 30 1-2. As poi nted out by e.g. lWalter et al.ld2003l 
2006) and llbarra et alj d2007i) . the non-existence of the Compton 
shoulder could be due to a strongly inhomogeneous absorbing 
medium. Since the strength of the shoulder is strongly depen- 
dent on the assumed accretion geometry, further work using self- 
consistent modeling of the absorption, fluorescent line formation 
and Compton shoulder formation is required. We will present 
such self-consistent analyses, as well as a more detailed study of 
the variability of the source, in a future publication. 
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